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We report electronic and magnetic structure of arsenic-doped manganese clusters from density- 
functional theory using generalized gradient approximation for the exchange-correlation energy. We 
find that arsenic stabilizes manganese clusters, though the ferromagnetic coupling between Mn atoms 
are found only in Mn2As and MruAs clusters with magnetic moments 9 fiB and 17 /is, respectively. 
For all other sizes, x = 3, 5-10, Mn^As clusters show ferrimagnetic coupling. It is suggested that, if 
grown during the low temperature MBE, the giant magnetic moments due to ferromagnetic coupling 
in Mn2As and Mn4As clusters could play a role on the ferromagnetism and on the variation observed 
in the Curie temperature of Mn-doped III-V semiconductors. 

PACS numbers: 61.46.+W, 36.40.Cg, 75.50.Pp 



Manganese is a unique 3d transition metal element due 
to its unusual and fascinating electronic and magnetic be- 
havior as atom, cluster, crystal and as well as impurity. 
The filled 4s and half-filled 3c? shell and the large energy 
gap ~ 8 eV between them prevent significant s — d hy- 
bridization and therefore they do not bind strongly as 
Mn atoms begin to form cluster [J, Similarly, a-Mn, 
the most stable form of bulk Mn, has the least binding 
energy among all the 3d transition metal elements. The 
magnetic properties of Mn clusters are particularly inter- 
esting. An early electron spin resonance study on small 
Mn clusters in an inert matrix suggested ferromagnetic 
ordering with magnetic moment 5 \ib /atom Q . More re- 
cently, through Stern-Gerlach (SG) molecular beam ex- 
periment , Knickelbein found that manganese clusters 
in the size range Mns - Mngg display ferrimagnetic order- 
ing with a maximum magnetic moment ~ 1.72 /is/atom 
for Mni2, despite of the fact that no known bulk phase 
of Mn displays such ordering. Recent density functional 
theory (DFT) calculations [H Q confirm this ferrimag- 
netic coupling in the size range x >5 in Mn x clusters. 

Manganese doped semiconductors, such as (GaMn)As 
and (InMn)As have attracted considerable attention be- 
cause of their carrier induced ferromagnetism The 
Mn dopants in these III-V dilute magnetic semiconduc- 
tors serve the dual roles of provision of magnetic mo- 
ments and acceptor production. Earlier experimental re- 
sults on Gai_ x Mn x As indicate a nonmonotonic behav- 
ior of Curie temperature T c (x), first increases with the 
Mn concentration x, reaching a maximum of 110 K for 
x ~5% and then decreases with the further increase of x. 
However, recent experimental studies 0,0 , under care- 
fully controlled growth and annealing conditions, suggest 
that the metastable nature and high defect content of 
low temperature molecular beam epitaxy (MBE) grown 
Gai-^Mn^As may be playing an important role in de- 
termining the magnetic properties. In particular, under 
suitable conditions, careful annealing could lead to an 



enhancement in T c with increasing Mn content [|| Q. 
Chiba et al. Q have reported a maximum T c as high as 
160 K, for 7.4 % Mn concentration, in a layered struc- 
ture. However, the detailed microscopic nature and the 
role of the metastable defects in Gai-^Mn^As are not 
yet well understood both theoretically and experimen- 
tally and demands more investigation. Strong segrega- 
tion tendency of the doped transition metal atoms into 
the semiconductor host [£| makes Mn clustering, around 
As in Gai-^Mn^As samples, a very important issue. 

In this Letter, we address the possibility of such clus- 
tering of Mn around As, and if so, then what is the nature 
of Mn-Mn magnetic coupling? We have, indeed, found 
that the binding energy of Mn clusters are substantially 
enhanced by single As doping by having their hybridized 
s — d electrons bond with p electrons of As. This stabi- 
lization is accompanied by the ferromagnetic or ferrimag- 
netic coupling between Mn atoms, all of which results in 
the large cluster magnetic moments, which could play a 
crucial role in the observed ferromagnetism and deter- 
mining the Curie temperature T c , in Mn-doped GaAs 
and InAs. 

Calculations have been carried out using DFT, within 
the plane wave method We have employed 

projector augmented-wave method ^{j and Perdew- 
Burke-Ernzerhof exchange-correlation functional |lll | for 
the spin-polarized generalized gradient approximation 
(GGA), as implemented in the VASP package [l^. The 
wave functions are expanded in a plane wave basis set 
with the kinetic energy cutoff equal to 337.3 eV and the 
calculations were carried out at the T point. The 3d, 4s 
for Mn and 4s, 4p orbitals for As were treated as va- 
lence states. Symmetry unrestricted geometry optimiza- 
tions were performed using quasi Newtonian and con- 
jugate gradient methods until all the force components 
are less than 0.005 eV/A. Simple cubic supercells are 
used with neighboring clusters separated by at least 12A 
vacuum regions. Several initial structures were studied 
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to ensure that the globally optimized geometry does not 
correspond to the local minima, as well as, for all clusters, 
we have explicitly considered all possible spin multiplici- 
ties to determine the ground state magnetic moment. 

We begin our discussions with pure Mn x (x <10) clus- 
ters, some of which (x <7) have also been studied pre- 
viously within the all-electron calculation and different 
levels of GGA fj. We have found ferromagnetic cou- 
pling for Mn2, M113 and Mn4 with magnetic moment 5 
fiB /atom, the Hund's rule value for the free atom. Emer- 
gence of ferrimagnctic coupling starts form Mns. The 
total magnetic moments of Mn x clusters corresponding 
to the ground state geometries, for x = 5-10, are 3, 8, 5, 
8, 7 and 14 [Ib respectively (see Tabic QJ, which are in 
excellent agreement with the recent SG experiment 0. 
Because of the half-filled 3d and filled 4s shell and a sub- 
stantial gap between them, we found Mn2 is a weakly 
bound van der Walls (vdW) dimer with bond length 2.58 
A and binding energy 0.53 eV/atom, as it is evident from 
the low experimental value [13j. Antiferromagnetic Mn2 
dimer is 0.53 eV higher in energy, whereas for Mn3, the 
frustrated antiferromagnetic structure is almost degen- 
erate (0.05 cV higher in energy) with the ferromagnetic 
structure. As x increases, binding energy increases mono- 
tonically [FigEfa)], reaches a value 1.94 eV/atom, 66% 
of bulk Q-Mn, for Mnio- We note, generally, the Mn-Mn 
bonds between like spins are larger (2.30-2.85 A) than 
that of between opposite spins (2.20-2.75 A). Several iso- 
mers lying close in energy to the ground state of Mn x 
were found which will be discussed elsewhere Q, in de- 
tail. 

In the Fig.??, we depict the ground state geometries 
of Mn^As clusters with their spin ordering. The MnAs 
dimer has much higher binding energy, 1.12 eV/atom, 
and much shorter bond length, 2.21 A, than the Mn 2 
vdW dimer. We have repeated the calculations of MnAs 
dimer, including Mn 3p as valence state an d have 
obtained an optimized bond length 2.22 A and binding 
energy 1.08 eV/atom with same total magnetic moment. 
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FIG. 1: (a) Plot of binding energy/atom with x, for Mn^ and 
Mn^As. (b) Two energy gains, A 1 and A 2 , are plotted with 
x. These energy gains are defined as, A 1 = — [E(Mn x As) — 
E(Mn x ) - E(As)] and A 2 = — [£(M ni As) - ^(Miii-.As) 
- E(Mn)}. 



Binding energy increases substantially to 1.63 eV/atom 
for the isosceles triangular Mn 2 As. The Mn-Mn dis- 
tance, 2.59 A, in Mn 2 As is almost equal to the Mn-Mn 
distance in M112, 2.58 A, and Mn-As-Mn bond angle is 
found to be 68°. As more Mn atoms are added, the 
structures take three-dimensional shape, and the deter- 
mination of ground state become a delicate task, as more 
than one geometric and/or magnetic isomers, very close 
in energy to the ground state are found, which will be 
presented elsewhere in detail. As mentioned earlier, this 
was done by minimizing the total energy of each clus- 
ter with respect to all possible spin multiplicities for ev- 
ery geometric structure studied. The presence of an As 
atom makes M113AS cluster tetrahedral, however, M114AS 
is a Mn4 tetrahedron with As at a face cap. Mi^As 
is the only cluster whose geometry differs significantly 
from a pure Mng cluster. Mng is octahedral, whereas 
MngAs is a pentagonal bipyramid, where the As atom 
is trapped in the pentagonal ring. As x in Mn^As in- 
creases, binding energy increases very slowly from 1.63 
eV/atom for Mn 2 As and tends to saturate to a value 
2.21 eV/atom for MnioAs (Fig^a)). The shortest Mn- 
As bond length increases from 2.21 A for MnAs to 2.46 
A for MnioAs, which is 4% and 2% shorter than the 
Mn-As distance in a-MnAs and Gai-^Mn^ As, respec- 
tively, whereas the shortest Mn-Mn distance decreases 
from 2.59 A for Mn 2 As to 2.23 for MnmAs. Generally, 
we find the same trend as seen in the pure Mn^ clusters 
that the bonds between Mn atoms of opposite spin to be 
somewhat shorter (2.20-2.60 A) than the bonds between 
Mn atoms of like spin (2.50-2.90 A), whereas all Mn-As 
distances vary between 2.20-2.60 A. All the Mn-As-Mn 
bond angles in these clusters vary in between ~ 60-70°. 
All the clusters in the Fig.?? and their respective iso- 
mers are magnetically stable i.e. both the spin gaps are 
positive: the lowest unoccupied molecular orbital of the 
minority (majority) spin lies above the highest occupied 
molecular orbital of the majority (minority) spin. These 
two spin gaps, Si and 82 [l5|, for Mn 2 As (0.83 and 1.34 
eV) and M114AS (0.89 and 1.14 eV) are the highest among 
all clusters. As x increases, 5i and 82 decrease to a value 
0.47 and 0.35 eV, respectively, for Mn 10 As. 

The next important issue is to see whether these Mn 
clustering around single As arc at all energetically favor- 
able or not. To understand this point, we calculate two 
different energy gains, A 1 - the energy gain in adding 
an As atom to a Mn^ cluster and A 2 - the energy gain 
in adding a Mn atom to a Mni-^As cluster. Fig^a) 
shows that due to the lack of hybridization between the 
4s and the 3c? electrons binding energy of pure Mn x clus- 
ters are very small. However, as an As atom is attached, 
the 4s 2 electrons of Mn interact with the 4p 3 electrons of 
As, which results in the substantial enhancement in the 
bonding (Fig^a)), and consequently, A 1 increases with 
x, which finally tends to saturate (Fig^b)). A 2 gives 
the number that how many Mn atoms can be bonded to 
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TABLE I: Total cluster magnetic moments fj, x of pure Mn^ 
and Mn x As clusters, corresponding to the ground state, for 
cluster size x < 10. 
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a single As atom, which is still significant, 2.65 cV, for 
MnioAs. These behaviors of A 1 and A 2 indicate that the 
Mn clusters around As are energetically favorable and we, 
therefore, argue that they are, likely to be, present in the 
low temperature MBE grown (GaMn)As/(InMn)As. 

Now we ask the next and most important question that 
what is the Mn-Mn magnetic coupling in these Mn^As 
clusters? The total magnetic moments of Mn x As clus- 
ters corresponding to the ground state geometries are 
given m Table |I| These large magnetic moments gen- 
erally arise from the ferrimagnetic coupling between the 
moments at Mn sites with the exceptions for Mn2As and 
MiuAs, where the magnetic coupling is ferromagnetic. 
We note, generally, no change in the nature of magnetic 
coupling between the Mn sites in Mn x As clusters from 
their respective pure Mn x . However, only for MnaAs, 
the Mn-Mn coupling behavior changes form ferromag- 
netic to ferrimagnetic, due to As doping, where a ferro- 
magnetic tetrahedral structure with total moment 12 \xb 
is found to be 0.12 eV higher. The nature of the mag- 
netic coupling is clear from their respective constant spin 
density surfaces, which are plotted for M112AS, M^As, 
M114AS and MnsAs in the Fig.??. The ground state mag- 
netic moments of Mn x As, for x=l, 2, 3, 5, 7, 8 and 10 
can be represented as (/i x — 1) fiB, whereas for x=A and 
9, it can be expressed as ([i x — 3) hb, where \i x is the 
total magnetic moment of the Mn x cluster correspond- 
ing to the ground state or the first isomer Q. Local 
magnetic moment, M, at each site can be calculated by 

AI = Ia iP^( r ) ~ Pl( r )]dr, where p T (r) and p ; (r) are up- 
spin and down-spin charge densities, respectively, and R 
[l6j is the radius of a sphere centring the atom. For MnAs 
dimer, the magnetic moment at Mn site, Mmu, and at 
As site, Ma s , are 3.72 fis and —0.26 fiB, respectively. 
This large negative polarization of the anion, As, is due 
to the strong p — d interaction. In the Mn2As cluster, Mn 
atoms are ferromagnctically coupled with Ma/„=3.79 [is 
each, whereas Ma s is —0.14 hb- Mn atoms in the M113AS 
take frustrated antiferromagnetic structure with Mmu = 
3.1, 3.1 and —3.9 hb and Ma s has a value —0.21 [1b- 
For Mn4As, Mn atoms are ferromagnetically arranged 
(Fig.??) with an average Af^/„=3.66 hb and are coupled 
antifcrromagnctically with the As atom, Ma s =— 0-22 hb- 
For MnsAs, Mun varies between 3.04-3.72 fis with a 



local magnetic moment —0.23 [is at As site. Polarized 
neutron diffraction study found a local magnetic moment 
of -0.23 ± 0.05 [i B at the As sites for NiAs-type MnAs 
[l7j . which is very close to the present values for MnAs 
- MnsAs clusters. These Mn-Mn magnetic behavior are 
unlike the previous study of nitrogen-doped Mn x clus- 
ters, a;=l-5, by Rao and Jena 0|, where Mn atoms were 
found to be coupled ferromagnetically for all sizes. We 
observe, the negative polarization of As, Ma s , decreases 
sharply to —0.08 [Ib for Mi^As and further decreases 
monotonically to — 0.02 /is for Mn§As, however it be- 
comes positive, 0.04 and 0.02, for MngAs and MnioAs, 
respectively. For all Mn x As clusters, x =6-10, Mn atoms 
are coupled ferrimagnetically, and Mmu of highly coor- 
dinated atoms are very small than those of the surface 
atoms and it varies between 0.8 - 3.7 \ib- 

In conclusion, though the As atom induces clustering of 
Mn, the individual magnetic moments of Mn couple fer- 
romagnctically only for M^As and Mn4As clusters and, 
however, ferrimagnetically aligned for all other sizes. Not 
only that, the ferromagnetic M^As and Mn4As are even 
more stable than other sized ferrimagnetic clusters, as 
they have largest spin gaps. We believe that Mn clus- 
tering, during the low temperature MBE growth, could 
be responsible for the ferromagnetism and reported high 
T c in (GaMn)As and should be taken into account to 
formulate an adequate theory of ferromagnetism in III- 
V semiconductors If Mn-doped GaAs, under care- 
fully controlled growth and annealing condition, contains 
Mn2As and/or Mn4As clusters, which have large cluster 
magnetic moments, would enhance the local magnetic 
moment and consequently could enhance T c . On the 
other hand, large clustering would yield low T c due to 
their small local magnetic moment. This study also sug- 
gest that the similar mechanism could enhance the T c 
of (InMn)As and further investigation of Mn-clustering 
around oxygen could explain the ferromagnetism and ob- 
served high T c of Mn-doped ZnO in the metastable phase 
|20j . EXAFS studies would be very useful to see whether 
these clusters of Mn around As are, indeed, present in 
these samples and the gas phase experiments involving 
Mn clustering in a As-seeded chamber can yield direct 
information on the magnetic behavior of Mn x As clusters. 
We hope, our study will encourage such experiments. 
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Figure Captions 

1. Ground state geometries and the corresponding 
spin ordering of Mn^As clusters, £=1-10. Blue and 
red balls represent Mnj and Miij atoms, respectively. 
Green ball represents the As atom. The bond lengths 
are given in A. Magnetic polarization of As, is negative 
for MnAs - MngAs and, whereas, positive for MngAs 
and MnioAs. Note, Mn-Mn interaction is ferromagnetic 
only for M112AS and Mn4As clusters and ferrimagnetic 
for all other sizes. 

2. (a) Plot of binding energy/atom with x, for Mn x 
and Mn^As. (b) Two energy gains, A 1 and A 2 , are 
plotted with x. These energy gains are defined as, A 1 
= - [^(Mn^As) - £(Mn x j - E(As)] and A 2 = - 
[E(Mn x As) - ^(Mni^As) - E(Mn)]. 

3. Constant spin density surfaces for M^As, M^As, 
Mn 4 As and Mn 8 As corresponding to 0.04, 0.04, 0.04 and 
0.02 e/A 3 , respectively. Red and blue surfaces represent 
positive and negative spin densities, respectively. Green 
ball is the As atom, which has negative polarization in all 
these structures. Note ferromagnetic (M^As and Mn 4 As 
(left panel)) and ferrimagnetic (M^As and MngAs (right 
panel)) coupling between Mn atoms. 
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